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Abstract 

The adhesion of flowing particles and biological cells over fixed collecting surfaces is vitally important in 
diverse situations and potentially controlled by small-scale surface heterogeneity on the particle. Differences in 
the behavior of patchy particles (flowing over uniform collectors) relative to the reverse case of uniform particles 
(flowing over patchy collectors) are quantified. Because a particle rotates more slowly than it translates in the 
shear field near a collecting surface, the effective interaction time of a patch on a particle is larger than that 
of a patch on the collector, suggesting distinct particle capture tendencies in each case. This paper presents a 
new computational approach to simulate the near-surface motion (rotation and translation) of particles having 
nanoscale surface heterogeneities flowing over uniform collectors. Small amounts of 10 nm cationic patches 
randomly distributed on a net-negative particle surface produced spatially varying DLVO interactions that 
were computed via the Grid Surface Integration (GSI) technique and then combined with hydrodynamic forces 
in a mobility tensor formulation. Statistical analysis of simulated trajectories revealed fewer extrema in the 
fluctuating particle-collector separation of heterogeneous particles, compared with the reverse system geometry 
of uniform particles flowing past a heterogeneous fixed surface. Additionally, the patchy particles were captured 
to a lesser extent on uniform surfaces compared with the case of uniform particles flowing above patchy collectors. 
Such behavior was dependent on ionic strength, with the greatest differences obtained near a Debye length of 
1 — 4 nm for the 2 a = 500 nm simulated particles. 


1 Introduction 

The adhesive capture of flowing particles on fixed surfaces is a critical step in the transport of bacteria in ground- 
water PH [331, in the initial stages of innate immune response [20], and in the manipulation of particles for 
applications ranging from mining to microfluidic sensors [28j. In such processes, flowing particles (or cells) are 
typically understood in terms of uniform properties hues ni]: surface charge, hydrophobicity/hydrophilicity bal¬ 
ance, and van der Waals interactions. Spatial distribution of surface functionality, “surface heterogeneity”, on 
these naturally-occurring and technologically relevant particles is consistently neglected, at length scales that vary 
from a few nanometers to the particle radius [38[ [32]. 

Patterned colloids such as Janus particles @5111] and colloidal molecules m have attracted intense interest for 
their potential to assemble into ordered periodic structures, not only for photonic applications IHEZi, but also, 
more fundamentally, to mimic the phase transitions of molecules [35[ HI 02]. The ability of engineered chemical and 
topographical patterns on these particles to direct their interactions suggest mm that the randomly distributed 
heterogeneities on naturally-occurring particles will influence their adhesion as well, even though the length scales 
of such particle features fall well below the colloidal dimensions [27[ 1 36] . 

The literature of the past two decades contain numerous efforts that model the adhesion of flowing particles 
on surfaces regularly-patterned with features that approach the lengthscales of the flowing particles m- The ad¬ 
hesion of flowing micron-sized particles on nano-scale heterogeneous collectors, however, has only been examined 
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recently mmmm- Surface heterogeneity, defined as a random arrangement of surface features, was shown to 
influence particle adhesion significantly, and it was demonstrated, for instance, that randomly-positioned clusters 
of nanoscale cationic patches on the collector can drive colloidal particle deposition onto net negative collecting sur¬ 
faces [27]. Even when each nanoscale cluster or “patch” of positive charge, positioned on the negative background 
surface, is individually too weakly attractive towards uniform flowing anionic particles, the concerted action of 
many of such patches can capture particles. Modeling methods that describe particle-collector interactions have 
evolved considerably over the past decade, from analytical [33 and numerical [a m computations of systems 
characterized by surface roughness, to the development of the surface-element integration (SEI) |6j and the more 
recent grid-surface integration (GSI) [2 Q7L [OH [12] techniques. While the SEI technique involves the discretization 
of only one interacting surface, in the now-accepted mum GSI technique surface forces are integrated over 
differential elements on two opposing objects of arbitrary shapes. 

In many recent computational studies, patchy particles are modeled as units composed of distinct “atoms.” 
Specific attributes are assigned to each atom, depending on whether it belongs to the ‘patch’ or ‘core’ surface area 
of the particle. Zhang and Glotzer|36] performed Brownian dynamic simulations to model the self assembly of 
particles patterned with patches positioned at specific locations. Polyelectrolytes of varying shapes are also shown 
by molecular dynamic computations to assemble into charged, patchy colloids m- 

While theoretical work on patchy particle interactions at the molecular scale abound, only a handful of stud¬ 
ies describe colloidal interactions of heterogeneous particles flowing near a planar wall and the effects of such 
heterogeneities on particle deposition. In contrast to the studies on heterogeneous fixed collectors, the system 
geometry comprising heterogeneous particles presents the additional complexity that, in shear flow, the particle 
rotates, exposing different areal elements to the collecting surface as it flows past. This rotation and turnover 
in the dominant surface region poses a computational challenge in the calculation of time-dependent colloidal 
interactions. Only a limited number of studies have addressed computationally related problems: Sphere-plate 
and sphere-sphere DLVO interactions for spheres patterned with roughness (topographical heterogeneity) were 
computed numerically [ZU39] and analytically m and found to agree with experimental measurements m • More 
recently, Chatterjee et al. focused on the deposition of micro-particles onto Janus and patchy spherical collec¬ 
tors [10] . Chemical heterogeneity was modeled by patterning the spherical collector with adhesion-favorable and 
adhesion-unfavorable surface properties assigned to the Janus collectors or with relatively wide alternating stripes 
(much wider than the particle radius) that covered the spherical surface. Deposition of the colloidal particles onto 
spherical heterogeneous collectors was analyzed as a function of collector orientation, collector patterning, and 
colloidal particle velocity. 

Similarly, the adhesion of leukocytes on the vascular endothelium was studied computationally by Korn and 
Schwarz [26], who modeled flowing leukocytes that were captured by the binding of randomly-positioned immobi¬ 
lized receptors. The adhesive receptors on the leukocytes, modeled as hard spheres, were similar to the randomly 
distributed nano-patches on flowing micro-particles. The leukocyte capture, however, is fundamentally different 
from the heterogeneous particle adhesion, since the receptors underwent chemical binding according to a reaction 
probability when in close distance to a surface ligand. Interesting motions of leukocytes, that ranged from free 
flow to firm adhesion, were presented as a function of the attachment and detachment rates of the cells to the wall 
in adhesion regime diagrams. The treatment lacked the spatially varying longer range colloidal forces of interest 
in the current study, but was, nonetheless, a break-through in its attention to patterned particle rotation in flow. 
As a result of this advance, the computations predicted interesting motion signatures such as particle rolling and 
skipping. 

The present work describes a new computational approach to calculate particle-surface interactions, trajectories 
and adhesion probabilities of flowing heterogeneous particles on fixed planar uniform collectors. The method 
expands the implementations of the Grid-Surface Integration (GSI) technique by discretizing the particle surface, 
in order to include on its surface a random distribution of patches. DLVO (electrostatic and van der Waals) 
interactions are computed within the GSI technique on the basis of parallel flat plate expressions pH EB, and 
flow-driven motion is obtained from the mobility matrix of the hydrodynamics problem m • Two heterogeneous 
particle-collector system geometries are, moreover, compared in detail. While one of such systems comprises 
electrostatically heterogeneous particles flowing past a uniform anionic collector, the alternative system consists of 
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uniform anionic particles flowing past an electrostatically heterogeneous collector, such that surface compositions 
in each system are reversed. The two configurations differ in their number of particle trajectory fluctuations 
and their collection probability due to the particles’ rotational motion, given by the rotational velocity times the 
particle radius, that is slower than its translational motion, determined by the translational velocity. It is found 
that patchy particles present less fluctuating trajectories and smaller collection probabilities than those of uniform 
particles flowing past heterogeneous collectors. The findings are presented in terms of statistical parameters of the 
particle trajectories and adhesion regime diagrams. 

2 Sphere Discretization 

In previous work, the collector surfaces were electrostatically and topographically heterogeneous, while the flowing 
colloidal particles were smooth, uniformly charged spheres. A discretization of the spherical surface was therefore 
not needed, and only the heterogeneous collector was partitioned into small areal elements. The accurate mod¬ 
eling of systems of heterogeneous particles, however, does require the discretization of the spherical surfaces into 
differential elements, each of which can be assigned distinct surface properties. 

In the results presented in this work, spherical surfaces are discretized into regions of equal area using Leopardi’s 
[ 129 ] recursive zonal EQual area Sphere Partitioning (EQSP) algorithm. A schematic diagram of a sphere discretized 
into IVp = 500 regions (elements) is presented in Fig. Qa). 

The accuracy of the EQSP algorithm is first verified by discretizing a 2a = 1 fj,m diameter sphere into N p = 
25963 elements, such that the individual element area is d A{ ~ 121 nm 2 , corresponding to an 11 nm x 11 nm 
patch for consistency with previous studies |5] motivated by experiments [27, [36]. For each surface element on the 
sphere, its projected area on an horizontal plane is computed from 

dSi = |n • ej_| dAi , (1) 

where dAi is the area of the sphere element i, dSi is its projected area on the plane, n is the vector normal to the 
sphere element and ej_ denotes the direction normal to the projection plane. The sum of the elemental projected 
areas dSi is within 0.0038% of the exact value of 2 S = 2ir a? (for two hemispheres). 

3 The GSI technique 

In the GSI technique, the total force or energy of interaction is obtained from a pairwise summation of interactions 
between differential areal elements on the collector and particle surfaces nan: 

P(h) ei • e z dS, (2) 

particle wall 

where P(h ) is the force or energy of interaction per unit area between a particle’s areal element (dA) and a 
corresponding element on the substrate ( dS ). The unit vector e z specifies the direction normal to the surface, and 
the unit vector e\ indicates the direction between areal elements on the particle and the substrate (see Fig{l|b)). 
The summations in Eq. ([2]) are performed over all of the discretized elements on the particle and collector. The 
colloidal energy of interaction per unit area ( P{h ) = U A (h)) between each pair of areal elements is obtained by 
summing van der Waals (vdW) and electrostatic double layer (EDL) components, each of them computed from 
analytical expressions nasi] derived for parallel plates. The corresponding forces per unit area ( P(h ) = F A {h)) 
are obtained by differentiation of the energies with respect to the separation distance and included in the mobility 
matrix formulation to yield particle trajectories Bain El- 

Within the GSI technique, each areal element on each interacting surface can be assigned distinct properties, 
such that GSI computations can be implemented for particles and collectors of arbitrary shape, topography and 
chemical properties. 
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Figure 1: (a) Illustration of a sphere discretized into N p = 500 equal area regions, (b) Schematic diagram of a 
uniform particle of radius a interacting with a uniform flat collector. The local and minimum particle-collector 
separation distances are denoted by h and D, respectively. Differential surface elements dA on the particle and dS 
on the collector are also indicated, as well as the unit vector ei that points in the direction of the collector surface 
element. The rotational and translational velocities are Q x and V y , respectively, and the flow shear rate is 7 . 
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Figure 2: Energy-distance profile for a uniform particle interacting with a uniform collector. The particle’s surface 
is either discretized into equal area elements (with the EQSP algorithm) or treated as a non-discretized uniform 
surface. In both cases, DLVO interactions are computed with the GSI technique. The simulation parameters are: 
2a = 1 n m, 7 = 25 s -1 , k -1 = 4 nm, A H = 5xlCT 21 , T co n ec tor = ^sphere = -25 mV. 

4 Uniform particle-collector systems 

A schematic diagram of a uniform particle-collector system, in which all interacting surface elements are assigned 
the same properties, is presented in Fig. 0b). Negatively charged, flowing particles of radius a interact with 
the flat uniform collector substrate, which is also negatively charged. The local and minimum particle-collector 
separation distances are denoted by h and D, respectively. Due to the linear shear flow with shear rate 7 , the 
particle translates in the y-direction with velocity V y and rotates around an axis parallel to the collector surface 
with a rotational velocity Q x . The origin of the Cartesian coordinate system is the left edge of the collector. For 
interactions of 2 a = 1 //m diameter particles, the collector length is L = 30/rm, while shorter collectors of length 
L = 20/im were simulated for interactions of 2 a = 500 nm diameter particles. A spherical coordinate system with 
origin in the sphere’s center, is also defined. The inclination angle is 9, 0 < 6 < ir, and the azimuthal angle is (f>, 
0 < cf) < 27r. The radial direction is normal to the particle surface. 

Electrostatic and van der Waals interactions for the uniform system illustrated in Fig. [l](b) are obtained by 
implementing the GSI technique, and including the discretization of either one or both interacting surfaces. For 
the case of a uniform spherical particle, all the EQSP-generated differential elements are assigned the same surface 
properties, which also equal those of the uniformly patterned collector. 

In Fig. [2| the energy profile of a smooth and uniformly charged 2a = 1 ym diameter particle interacting with a 
flat surface is shown. Colloidal interactions are characterized by the Hamaker constant Ah = 5xl0~ 21 J, which is 
fixed throughout this study, and the inverse Debye screening length k^ 1 = 4 nm. The surface loading 0, defined as 
the area of the collector containing attractive (cationic) functionality relative to the total collector area, is equal to 
zero for both interacting surfaces, which carry an electrostatic potential T = —25 mV. The energy profile is thus 
defined by the repulsive electrostatic double layer (EDL) interactions between the uniformly and equally charged 
particle and collector surfaces. 

The interactions presented in Fig. [2] are computed by implementing the GSI technique and either including 
(solid line) or not (dotted line) the EQSP algorithm for the modeling of the uniformly charged, smooth sphere. 
As in Sec. [2j the uniform 2a = 1 yrn diameter particle is discretized into N p = 25963 areal elements, such that 
the area of each element is dAj ~ 121 nm 2 . In the case for which the spherical surface is not discretized into 
differential elements, the collector surface is partitioned instead. The collector grid consists of 91 square elements 
that represent a length of 2a = 1 yrn, such that the length of each square element is ~ 11 nm. The spherical 
and planar discretization schemes are specifically chosen so as to define equal-area elements, in this case of dAj ~ 
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Figure 3: Schematic diagrams of particle-collector systems with one heterogeneous surface, (a) Patchy sphere- 
uniform collector, (b) Uniform sphere-patchy collector. 


121 nm 2 . In heterogeneous systems, the size of one areal element can sometimes equal the surface feature size 
[5j though it should be smaller if the surface chemical patterns become relatively large. Discretization schemes, 
therefore, vary with particle and surface features dimensions. 

It is seen in Fig. [2] that both computational techniques yield results that are in perfect agreement, thus 
validating the implementation of the GSI-EQSP technique. 

5 Random distribution of cationic patches on the sphere. 

Particles and collectors with randomly positioned attractive cationic patches can be constructed by numbering 
the discretized sphere/collector elements and then choosing the location of each patch sequentially by mapping a 
uniform deviate. In this section, interactions between patchy particles and uniformly charged, flat collectors are 
studied for different particle sizes, surface loadings and Debye lengths. Detailed comparisons to interactions in 
uniform particle-patchy collector systems are also presented. The patchy particle and patchy collector systems are 
schematically depicted in Figj3j Patches on either heterogeneous surface are assigned a potential of 'I'het = 50 mV, 
while the other regions of the heterogeneous surface and the homogeneous surface bear a uniform electrostatic 
potential of ’Puni = —25 mV. The interactions between a patch and the uniform surface are therefore attractive, 
while the interactions between other regions are repulsive. The individual patch area equals that of the surface 
element (A patc h ~ 121 nm 2 ), such that the patches model 11 nm squares. Due to the Poisson distribution used to 
locate patches on both heterogeneous surfaces (planar or spherical), there are no regions on either heterogeneous 
surface where patches are preferably assigned. 

For a fixed area fraction of cationic patches of 0 = 0.17, trajectories of patchy particles flowing over uniformly 
charged collectors and those of uniformly charged particles flowing over patchy collectors are presented in Fig. 
|4](a)-(b) for two particle sizes and Debye lengths. For each case, the ratio of the particle’s rotational to translational 
velocities are plotted as a function of the horizontal displacement in Fig. m -(d). For clarity, Figj4](c) and the 
inset in such figure are presented separately in the Supporting Information. 

In particle trajectory and velocity calculations presented in Fig. [4j the particle is assumed to contact the 
collector at an arbitrarily small separation distance D = 5 = 1 nm that models the surface roughness Haim. 
When the separation distance D < 5, the particle either rolls in contact with the collector, or, arrests due to 
friction forces. The condition for particle arrest is found from a force balance in the direction of flow and reads 

m 

F y + — F r f < 0 (3) 
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Figure 4: Trajectories and velocity ratios for patchy and uniform particles of sizes 2a = 0.5, l^rm interacting at 
Debye lengths k -1 = 1,5 nrn, for a fixed surface loading 0 = 0.17. Rolling friction is computed with a friction 
coefficient of /xr = 1.3 x 10 - 4 m- (a)(b) Particle trajectories of patchy and uniform particles, for particle sizes 
of 2 a = 1 /jrn (a) and 2 a = 500 nm (b). (c)(d) Angular to translational velocity ratios of patchy and uniform 
particles, for particle sizes of 2 a = 1 /iin (c) and 2a = 500 nm (d). For clarity, all solid lines represent the patchy 
particle system, while all dotted lines indicate results obtained for the uniform particle (patchy collector) system. 
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where F y and T x are the shear-induced force and torque mm, respectively, and F r f = HrF^ is the rolling 
resistance, determined by a coefficient of rolling friction /ur and the normal force F^ = —Frlvo- The particle 
therefore arrests on the collector once the DLVO forces, multiplied by an appropriate friction coefficient, are 
sufficiently attractive to overcome the shear-induced effects, or, when 


F 


y H-- + VrFdlvo < 0 • 


(4) 


The rolling friction is computed with a friction coefficient of /j,r = 1.3 X 10 - 4 eh. 

In Fig. 0a), the upper lines represent trajectories of 2 a = 1/rrn diameter patchy and uniform spheres translating 
in shear flow above uniform and patchy collectors respectively, for a Debye length of k _1 = 5 nm. Due to the large 
Debye length, both particles translate at a relatively large separation distance from the collector, and the particles 
do not contact or deposit on the collector. The average separation (mean ± standard deviation) for the patchy 
particle is D/a = 6.64 x 10~ 2 ± 2.08 X 1CV 3 , while that of the uniform particle is D/a = 6.60 x 10 -2 ± 2.13 x 10 -3 , 
in agreement with the secondary minimum in the energy profile D m ; n /a = 6.62 x 10 -2 , computed with the GSIua 
technique [5] . Although the mean separation distance and its standard deviation are consistent for the patchy 
and uniform particles, collection probabilities are found to be different in patchy particle systems than in patchy 
collector systems. This behavior is examined in more detail in Secs. 5.1 and 5.3 In Fig. 0c), the ratio of the 
rotational and translational velocities for a Debye length k^ 1 = 5 nm are shown by the lower solid and dotted 
lines, which indicate results for the patchy and uniform spheres, respectively. In both cases, the ratios fluctuate 
around the same values, that are less than unity when the particle is not in contact with the collector, confirming 
that the particles translate faster than they rotate. 

The lower solid and dotted lines in Fig. 0a) show trajectories of 2 a = 1 /jrn diameter patchy and uniform 
particles, respectively, for a Debye length = 1 nm. It is seen that both particles frequently contact the 
collector surface (at an arbitrarily small distance <5 = 1 nm taken as representative of surface roughness |14] 1 to 
yield a trajectory that is characterized by alternating periods of rolling motion and of free flow in close proximity 
to the collector. Both particles maintain a separation distance D < 5 nm as they flow above the entire simulated 
collector. Velocity ratios as a function of the horizontal displacement are shown by the upper solid and dotted lines 
(for patchy and uniform particles respectively) in Fig. 0c). The rolling periods for both particles are identified by 
velocity ratios that are equal to unity, as it is assumed that the particle does not slip when it contacts the surface 
m , while the translation between those periods corresponds to velocity ratios that fall below unity. 

Fig. @b) presents trajectories of 2a = 500 nm particles at Debye lengths k 1 = 1,5 nm. The patchy spheres 
are discretized into N p = 6489 elements, such that the surface area of each element remains d A{ ~ 121 nm 2 , 
and the patchy collectors are modeled with a grid consisting of 45 square elements that represent a length of 
2a = 500 nm, such that the length of each square element is ~ 11 nm. At a large Debye length of kT 1 = 5 nm, the 
patchy particle (upper solid line) translates at a separation D/a of (0.134 ± 5.8 x 10 -3 ), while that of the uniform 
particle (upper dotted line) is (0.133 ±6.72 x 1CP 3 ), in good agreement with the secondary minimum in the energy 
profile of D m i n /a = 0.134, computed with the GSIua technique [5]. The amount of spatial fluctuations of the 
uniform particle trajectory is larger, as expected from the larger standard deviation of the separation distance. 
More importantly, the uniform particle (upper dotted line) adheres on the patchy collector, while the patchy 
particle (solid upper line) flows above the entire uniform collector, without adhering. The respective velocity 
ratios, presented in the lower solid and dotted lines in Fig. 0d) are smaller than 1. Both particles translate fast, 
with respect to their rotational motion, and, for the uniform particle, the ratio increases to 1 when the particle 
is arrested on the collector due to a rolling resistance in the direction of flow. Such rolling resistance retards the 
particle’s motion as a consequence of elastic deformations of the surfaces in contact HU. 

For a Debye length of k 1 = 1 nm, trajectories of patchy and uniform 2 a = 500 nm diameter particles are 
shown by the lower solid and dotted lines in Fig. EDM, and enlarged for clarity in the inset of the same figure. Both 
particle trajectories exhibit rolling periods that alternate with free flow segments in which the particles translate 
in close proximity to the collector. Due to friction forces, both particles ultimately adhere on the collector when 
the DLVO attraction becomes sufficiently large. The velocity profiles, denoted by the upper solid and dotted lines 
in Fig. 0d), reach a value of unity, thus indicating rolling periods, but fluctuate around values lower than unity 
in the trajectory segments in which the particles loose contact with the collector. 





Figure 5: Collection probability curves of 2 a = 500 nm diameter, patchy particles adhering to a uniform anionic 
collector and of uniform anionic particles adhering to a patchy collector. 0 is the area fraction of positive patches 
on either patchy surface. 


5.1 Collection Probability 

In Sec. [5j a few characteristic examples of particle trajectories were computed for different patchy/uniform particle 
sizes and Debye lengths. The results presented in this section are statistically significant collection probabilities 
obtained from a large number of computed trajectories. 

In the present study, collection probabilities are defined as the ratio of adhered particles (successes) to the 
total number of simulated particle trajectories (trials), 


V = 


N d 

IVtotai 


(5) 


where N d is the number of particles deposited on the collector and IV tota i is the total number of particle trajectories 
simulated. The collection probability estimate given by Eq. ([ 5 ]) is the ratio of successful adhesion attempts to 
the total number of attempts and resembles adhesion probabilities |25] or the available surface function (ASF) 
defined as the normalized particle adsorption probability in the limit N at tempts —> 00 [3]- It should be noted that, 
in the present study, particle trajectories were computed sequentially for collectors with no adhered particles and 
particle adhesion is thus determined by particle translation and rotation in shear flow subject to the computed 
electrostatic and van der Waals interactions. 

In Fig. [5j collection probability curves are presented for patchy and uniform particles at a Debye length 
At -1 = 4 nm. The total number of particle trajectories simulated for each point on this plot is IVtotai = 1000. 
Wilson score intervals [53] are constructed to provide more statistically relevant predictions and the error bars 
in Fig. [ 5 ] correspond to the 95% confidence interval. Each of the trajectories computed for patchy particles was 
performed for a different heterogeneous particle in which the patches are randomly distributed on the spherical 
surface. Alternatively, uniformly charged particles were simulated to flow over distinct heterogeneous collectors, 
each of which was patterned with patches located at randomly chosen locations. Once the sphere or planar surface 
was discretized into elements with equal area, the algorithm used to choose the location of a patch was identical. 

Unlike deposition rate curves, the collection probability curves presented in Fig. [5] do not contain any rate 
information. In similar computational studies performed on the basis of the Random Sequential Adsorption (RSA) 
model, the particle adsorption kinetics is characterized by plots of the absorbed particle coverage as a function 
of a dimensionless adsorption time that depends on the number of simulated adhesion attempts [a e], which 
is most relevant when a significant fraction of the surface is covered by adhered particles. Particle deposition 
modeled with RSA simulations does not include the effects of DLVO interactions, but instead depend on a series 
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of attempts in which locations on the heterogeneous surface are randomly chosen. A particle will be considered to 
be irreversibly adhered at a specific location if an unoccupied adhesion site is present at such location, and if the 
additional particle would not overlap with previously adhered ones. Experimental work [22 El EH |36] involving 
heterogeneous collectors with nanoscale patches has been focused on adhesion thresholds and initial deposition 
rates (when the fraction of the collector covered by adhered particles is negligible). In the present study, each 
particle was thus sampled one at a time and each adhesion attempt is represented by a different (uniform or 
patchy) particle which adheres, or not, on the (patchy or uniform) collector, which is not occupied by previously 
adhered particles. Particle adhesion is then determined by DLVO particle-collector interactions, as the particle 
translates and rotates in shear flow. 

The collection probability curves presented in Fig. [5] provide an estimate of the particle adhesion threshold, 
defined as the smallest surface loading for which particles adhere on the collector (@ c = min(@) : rj > 0). 
Interestingly, for the chosen set of parameters, the adhesion threshold for the patchy sphere on a uniform collector 
(0 C = 0.14) is larger than that of the uniform sphere on a patchy collector (0 C = 0.12). The higher threshold 
for the patchy sphere is attributed to its decreased tendency to contact and adhere on the uniform collector, with 
respect to the uniform sphere flowing over the patchy collector. The patchy particle is less prone to adhere on the 


uniform collector due to its larger patch interaction time, as described in detail in Sec. 5.3 


5.2 Number of local extrema in particle-collector separation 

Trajectories for the patchy particle and patchy collector systems differ in the extent of fluctuations in the particle- 
collector separation. Such differences can be quantified by computing the number of local extrema in the trajecto¬ 
ries as a function of the surface loading 0, for each particle type (patchy or uniform). Due to the inherent random 
nature of the heterogeneity, the number of extrema for each particle type at each value of 0 is obtained as an 
average performed over many trajectories. 

Since not all particle trajectories denote particle adhesion, the average density of extrema for a fixed 0 is 
obtained as the linear combination 


< -Next >=p < iVext > a +(1 ~ p) < N ext >na 


( 6 ) 


where the subscripts a and na denote ‘adhesion’ and ‘no adhesion’ respectively, 


Np 

-N-total 


(7) 


Np is the number of particles deposited on the collector, and AI to tai — 200. Eq. ©> is derived in detail in the 
Supporting Information. The variance of < N ex t >, 


2 2 2 
^Next = P ^Ne 


+ (i -p)V 


Next, r 


( 8 ) 


is obtained from an error propagation expression derived for Eq. ([b]), assuming that, for the specific case of the 
trajectories considered in these results, p is a precisely known constant. 

In Fig. [6j the number of extrema is plotted as a function of the surface loading 0 for n~ l = 4 nm. It is 
seen that for both patchy and uniform particles the average number of local extrema remains relatively constant 
for low values of the surface loading of patches 0, but increases as the surface loading is further increased. As 
described in Sec. [5j interactions with patches are defined as attractive. For increasing surface loading of patches, 
thus, the increasingly stronger attractive interactions decrease the particle-collector separation distance, resulting 
in a larger number of spatial fluctuations (number of extrema) in the particle trajectories. The average number 
of extrema starts increasing at 0 ~ 0.17 and at 0 ~ 0.19 for the uniform and patchy particles, respectively. The 
standard deviation in the number of extrema increases as well, and more significantly than the number of extrema 
itself, with increasing surface loading 0. As the number of adhering particles increases, the existence of different 
adhesion sites on each heterogeneous collector multiplies the number of different trajectories that denote particle 
adhesion, thus increasing the standard deviation in < N ext > at larger values of 0. As shown in Fig. [6j the 
standard deviation increases dramatically when the adhesion probability approaches unity. 
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Figure 6: Weighted number of extrema as a function of © for patchy and uniform 2 a = 500 nm particles, averaged 
over Wot — 200 trajectories. 


The total number of extrema for the patchy particle is smaller (Wot,p — Wot,u/2) than that of the uniform 
particle. As described in detail in Sec. |5.3[ the larger patch interaction time of the patchy particle partially 
precludes interactions with other heterogeneity regions, which gives rise to smoother trajectories, or, equivalently, 
a smaller number of extrema. For each particle type (patchy and uniform), the maximum and minimum number of 
extrema are also computed and found to be essentially equal, A r maXi p (,,) = JV min iP ( u ) = Wotai, p(u)/2, as expected 
for a fluctuating trajectory line. 


5.3 Patch interaction times 


Differences in the trajectories and adhesive behaviors of patchy and uniform particles can be interpreted in terms 
of an effective time period over which one patch influences the electrostatic double layer interaction between a 
particle and the collector. 

For simplicity, the patch interaction time is defined as the time a sphere surface element interacts within the 
zone of influence (ZOI) on the collector. The radius of the ZOI is defined approximately by the overlap of the 
Debye layers [171136] and given by i?zoi = 2\/« 1 a. The newly defined patch interaction time is, in fact, the 
maximum patch interaction time, since the particle could adhere on the collector by translating a distance that is 
smaller than the Rzoi- For a uniform particle translating above a patchy collector, the patch interaction time is the 
time required for the particle to translate a distance of 2i?zoi- For a patchy particle, translating above a uniform 
collector, the patch interaction time is the time required for the particle to rotate through the corresponding 
subtended angle. 

A particle in shear flow in close proximity to a planar surface is shown schematically in Fig. [7](a) . The particle 
translates with velocity V y and rotates at an angular velocity fl x . For the case in which the patch is located on 
the collector, the appropriate interaction time depends on the translational velocity V y , as illustrated in Fig. 00. 
and is defined as 


Ttr = 


2-Rzoi 

Vv 


(9) 


If the patch is located on the spherical surface, as depicted in 
rotational velocity 

2 a 


Trot — 


D, 


Fig. [7](c), the interaction time is a function of the 

( 10 ) 


where 2 a is the angular displacement that corresponds to a linear displacement of 2i?zoi> and sin(a) = -Rzoi/(« + 
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(a) 



Figure 7: (a) Schematic definition of the radius of the Zone of Influence Rzoi- Debye layers of width k _ 1 around 
each interacting surface, particle velocities Da, and V y , and the angular displacement a that corresponds to a linear 
displacement of Rzop are also indicated. (b)(c) Schematic diagrams illustrating the appropriate displacements and 
velocities that define patch interaction times for patchy collectors and spherical particles, (b) Linear displacement 
and linear velocity V y for the case of a patchy collector (uniform particle), (c) Angular displacement and rotational 
velocity £l x for the case of a patchy particle. 


.-ii 


k '). The “rotational” interaction time defined by Eq. (10) is equivalent to a translational interaction time for 
small angular displacements (sin(a) ~ a) defined in terms of an effective translational velocity, V y = afl x , 


Trot 


2i?zoi ^ 2Rzoi 
(a + K -1 )!^ aQ x 


2i?ZOI 



( 11 ) 


For particle sizes and Debye lengths frequently chosen in experimental studies PE2U23], the angular displace¬ 
ments are indeed small. For example, the approximation sin (a) ~ a presents an error of less than 1% for particle 
sizes 2 a = 500 nm - 1 /mr and Debye lengths k - 1 = 2 - 5 nm. When the particle is not in contact with the 
collector, it is important to note that V y < V y , as the particle rotates more slowly than it translates H2HB1. The 
equality V y = V y holds, however, when the particle rolls without slipping m, as shown in Figs. |4j(c)-(d). 

Computations of patch interaction times are most relevant for particle trajectories at surface loadings for which 
the adhesion probability is small, i.e. near the adhesion threshold. The secondary minimum in the energy-distance 
profile is chosen as the appropriate separation distance for the computation of the velocities required in Eqs. (§- 
(10) and is computed using the GSIjja technique [5]. The rotational and translational velocities are obtained 
from the mobility matrix formulation of the hydrodynamics mmm and the hydrodynamic functions within 
such matrix are computed with the functional forms given by Duffadar and Davis m at a separation distance 
h = D. 


sec. mm. • 


In Fig. [8j patch interaction times are presented as a function of 0 for patchy and uniform particles of varying 
size and at Debye lengths k,~ 1 = 2,4 nm. In all cases, the patch interaction time slightly increases with 0 because 
stronger attractive interactions decrease the particle-collector separation, which reduces the particle velocities. 
The values of 0 considered in patch interaction time calculations are those for which a secondary minimum in the 
energy-distance profile exists. 

For a fixed particle size of 2 a = 500 nm, patch interaction times as a function of 0, for both particle types 
(patchy, uniform), and two Debye lengths k v = 2,4 nm are presented in Fig. |8|a). At a constant Debye length of 
k~ 1 = 4 nm, the interaction times of the patchy particles are significantly larger than those of the uniform spheres 
(in average, for the presented range of 0, r ro t — 1.75 Tt r ). The same trend is observed for the lower Debye length 
K -1 = 2 nm. The larger patch interaction time of the patchy particle limits heterogeneous (attractive) interactions, 
because patchy particles interact with each uniform collector element for longer periods of time. Uniform particles, 
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Figure 8: Dimensionless patch interaction times vs. surface loading, for patchy and uniform spheres of varying 
particle sizes and at different Debye lengths, (a) Patch interaction times of patchy and uniform particles at Debye 
lengths k -1 = 2,4 nm. The particle size is fixed at 2 a = 500 nm diameter, (b) Patch interaction times of patchy 
and uniform particles, of sizes 2a = 500 nm and 2 a = 1/xm. The Debye length is fixed at = 4 nm. (c) Ratios 
of translational to rotational patch interaction times for varying particle sizes and Debye lengths. 
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however, interact with each patch on the collector for shorter periods of time, allowing for an increased number 
of interactions with multiple patches within the same time period. In agreement with all the results presented 
within this work, the larger patch interaction time for patches on the particle translates to a smaller number of 
trajectory fluctuations (number of local extrema) and lower collection probabilities. It is thus also suggested that 
adhesion thresholds would be higher for patchy particles. 

For either placement of patches, on the particle or on the collector, it is also shown in Fig. |8](a) that, for a 
fixed particle size, patch interaction times are larger for the larger Debye length. At a constant particle size, a 
change in the Debye length has two competing effects. An increase in the Debye length corresponds to an increase 
of i?zoi = 2\/ n~ 1 a, suggesting the interaction time should increase at larger values of k - 1 . At the same time, 
however, a larger Debye length moves the secondary minimum farther away from the collector, which increases 
the particle velocity. At a constant particle size of 2 a = 500 nm, the increase in the size of the ZOI overcomes the 
increase in the particle velocity, such that the patch interaction times increase with increasing Debye length. 

The dependence of the patch interaction time on the surface loading 0 at a constant Debye length of k^ 1 = 4 
nm is shown in Fig. |8](b) for patchy and uniform particles of sizes 2 a = 500 nm and 1 fira. For both particle 
sizes, the patch interaction times of the patchy particle are larger than those of the uniform particles because the 
rotational velocity is significantly smaller than the translational velocity. The effect of particle size on the patch 
interaction time, at a constant Debye length, is dual, just as the effect of the Debye length at a constant particle 
size, previously described. At a fixed Debye length, the larger particle has a larger ZOI, thus, it could be expected 
that the larger particles will exhibit the larger interaction times. It is the smaller particle, however, for which the 
interaction times are larger. The smaller ZOI of the smaller particle defines more localized, attractive interactions 
that reduce the flowing particles’ velocities, to yield larger patch interaction times. 

The ratio of the translational and rotational patch interaction times for the particle sizes and Debye lengths 
presented in Figs. [8](a)-(b) are shown as a function of 0 in Fig. [8](c). It is readily noted that 


Ttr 

Trot 


2Rzoi/V y af2 a 


2R zol /aQ x V y 


( 12 ) 


The velocity-, or equivalently, patch interaction time-ratios given by Eq. (12) and shown in Fig. [8](c) are smaller 
than 1, and thus indicate that the particle translates faster than it rotates, in agreement with the results presented 
in Figs. |4](c)-(d) for trajectory segments in which the particle does not contact the collector. It is shown in the 
Supporting Information that the velocity ratio a^l x /V y increases as the separation (D/a) —» 0, such that smaller 
separations (due to stronger attractions at larger values of 0) lead to higher ratios of the interaction times, as 
shown in Fig. [8](c). 


5.4 Adhesion regime diagrams 

Interesting particle motions can be distinguished on the basis of system parameters and then presented in adhesion 
regime diagrams p6, H]. Duffadar and Davis m presented adhesion regime diagrams as functions of the Debye 
length and surface loading 0 in order to distinguish the ‘arrest’, ‘skipping and rolling’, and ‘no arrest’ behaviors 
of micro-particles in shear flow above a patchy collector. 

Similarly, an adhesion regime diagram is presented in Fig. [9j on the basis of collection probabilities, such as 
those computed in Sec. 5.1 plotted as a function of k - 1 and 0. The diagram presents adhesion/no adhesion 
regimes for 2 a = 500 nm patchy and uniform spheres. As expected from the results presented in previous sections, 
adhesion thresholds estimated from collection probability curves are consistently higher for patchy particles than 
for uniform ones, such that the adhesion regime of the patchy sphere is smaller than that of the uniform sphere. 

The differences in thresholds are more noticeable at relatively low Debye lengths. At a high Debye length of 
k~ 1 = 10 nm, computations within this work revealed an equal threshold of 0 C = 0.19 for both particle types. 
Due to a larger ZOI, the heterogeneity is smeared over a larger (planar or spherical) area and the interactions 
assume a more mean-field like character, irrespective of their specific location on the collector or sphere. 

It is interesting to note that each of the ‘adhesion’ and ‘no adhesion’ regimes can be subdivided to account 
for rolling/skipping motions and thus distinguish between 4 typical dynamic behaviors, given by the adhesion/no 
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Figure 9: Adhesion regime diagram, presented as the dependence of the adhesion thresholds 0 C estimated from 
collection probability curves on the Debye length for 2 a = 500 nm patchy and uniform particles. 


adhesion and surface contact/no surface contact possible combinations. Particles of 2 a = 500 nm diameter, with a 
Debye length = 5 nm (upper lines in Fig. Qb) are examples of particles that do not contact the collector and 
do not adhere on the surface (patchy sphere) and that do not contact the collector but do adhere on the surface 
(uniform sphere). Alternatively, at a low Debye length k - 1 = 1 nm, small patchy and uniform spheres contact 
the collector and adhere on it (lower lines in Fig. [4)(b) ), while large particles contact the collector but without 
adhering (lower lines in Fig. |4](a)). 

6 Conclusion 

Particle-collector systems with one heterogeneous surface, which is either the spherical particle or the planar 
collector, are characterized in detail for different system geometries and a range of particle sizes, Debye lengths and 
electrostatic potentials. DLVO interactions are computed by implementing the GSI technique, and incorporated 
in a mobility matrix formulation of the dynamics problem that yields particle trajectories as the particle translates 
in shear flow above a flat collector. 

The (particle)-(patchy collector) system, thoroughly studied in previous work, does not require spherical sur¬ 
faces to be discretized into areal surface elements. To model interactions of systems that include patchy particles, 
however, the inclusion of such a discretization scheme in the computational model is, in fact, essential. The use of 
the recursive zonal EQual area Sphere Partitioning (EQSP) algorithm for the discretization of spherical particles 
into small, equal-area elements is validated through energy computations of the uniform particle-collector system. 
DLVO interactions for (patchy particle)-(collector) systems are thus modeled by incorporating EQSP-generated 
spherical surface elements within the GSI technique. The system’s dynamic behavior is also obtained in this case 
from mobility matrix computations. 

Differences in the adhesive and dynamic behaviors of particle-collector patchy systems, in which only one 
surface is patterned with nanoscale features, are quantified by computations of collection probabilities and of 
average numbers of trajectory local extrema. Patch interaction times in each case are also defined. 

The lessened tendency of the patchy particle to adhere on a uniform collector, with respect to that of the 
uniform particle adhering on a patchy collector, is attributed to larger patch interaction times for the patchy 
particle. A larger interaction time precludes multiple interactions with many heterogeneous surface elements in a 
given time period, which translates into fewer attractive interactions. Moreover, larger interaction times reduce 
the amount of spatial fluctuations exhibited by patchy particles interacting with uniform collectors, and lead to 
a smaller number of local extrema in the trajectory, as shown in the results presented in this work. Spatial 
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variations in the trajectory, in turn, correlate with the extent of interactions with heterogeneous surface elements 
and ultimately provide insight into the adhesive character of the system for a given set of parameters. Adhesion 
thresholds estimated from collection probability curves of patchy particles are indeed larger than those of uniform 
particles adhering on patchy collectors. 

In summary, a new computational approach is introduced in this work to compute particle trajectories, along 
with statistical measures that characterize those trajectories, of patchy particles flowing above uniform planar 
collectors. The EQSP discretization scheme of spherical surfaces is incorporated in GSI computations of DLVO in¬ 
teractions and in the mobility matrix formulation that yields particle trajectories. The use of this newly developed 
simulation technique can be extended, for instance, to model particle-collector systems of two heterogeneous sur¬ 
faces, and of surfaces covered with many types of heterogeneity, such as patches, pillars, and spring-like structures 
that resemble polymer brushes or cellular receptors, which will be pursued in future work. 
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